function sim = set_up_mit(eq_SS, param, glob, options)
    
    
    sim.C = ones(options.T,1);
    sim.D = ones(options.T,1);

    sim.A          = zeros(options.T, 1);
    sim.A(1)       = -.0163; %
    sim.A(2:100)   = sim.A(1)*(.685.^(2:100));
    sim.A          = exp(sim.A);
           
    sim.M          = param.M*ones(options.T,1);

    
    sim.C = sim.A;
    
  %  sim.D          = eq_SS.D*ones(options.T,1);%.*((sim.cE./param.ce));
  %  sim.D(2)       = Din(2);
  %  sim.D(3:100)   = sim.D(end) + (sim.D(2) - sim.D(end))*.1.^(3:100);
    
    sim.W          = sim.C;
 
    sim.cE         = param.ce*ones(options.T,1);
 
    sim.cF         = param.mu_f*ones(options.T,1);
    
    sim.mu         = nan(size(eq_SS.L,1),options.T);
    sim.mu(:,1)    = eq_SS.L;
    
    sim.y          = nan(size(eq_SS.L,1),options.T);
    sim.p          = nan(size(eq_SS.L,1),options.T);
    sim.bp         = nan(size(eq_SS.L,1),options.T);
    
    sim.enter      = nan(size(glob.agrid,1), options.T);
    sim.exit       = nan(size(eq_SS.L,1),options.T);
    
end